Hrmed  Services  lechnical  Information  Hgency 

Because  of  our  limited  supply,  you  are  requested  to  return  this  copy  WHEN  IT  HAS  SERVED 
YOUR  PURPOSE  so  that  it  may  be  made  available  to  other  requesters.  Your  cooperation 
will  be  appreciated. 


NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


Reproduced  by 

DOCUMENT  SERVICE  CENTER 

KNOTT  BUILDING,  DAYTON,  2,  OHIO  ' 


I 


(ttotomhia  Hutoraitg 

in  tlyr  (Etlg  nf  5fan  fork 


go 


DEPARTMENT  OF  CIVIL  ENGINEERING 


FURTHER  STUDIES  OF  THE  RESPONSE  OF  A 
CYLINDRICAL  SHELL 
TO  A TRANSVERSE  SHOCK  WAVE 


M.  L.  BARON  and  H.  H.  BLEICH 


Office  of  Naval  Research 
Contract  Nonr-266(08) 
Technical  Report  No.  10 
CU- 1 0-  5 3-ONR-266 (08 ) -C.E. 
December  1953 


(Mumbta  IntwrHttg 

in  tlf?  (City  of  2mu  $orfe 


DEPARTMENT  OF  CIVIL  ENGINEERING 


FURTHER  STUDIES  OF  THE  RESPONSE  OF  A 
CYLINDRICAL  SHELL 
TO  A TRANSVERSE  SHOCK  WAVE 


by 


M.  L.  BARON  and  H.  H.  BLEICH 


Office  of  Naval  Research 
Contract  Nonr-266(08) 
Technical  Report  No.  10 
CU- 1 0-  5 3-ONR-266  ( 08 ) -C.  E. 
December  1953 


TABLE  OF  CONTENTS 


Page 

List  of  Symbols  iii 

I  Introduction v 

II  Interaction  between  Shell  and  Medium  1 

III  Response  of  Shell  to  a Step  Shock  Wave  for  n ^ 1 3 

IV  Accuracy  of  the  Approximation,  Rq.  (11-20)  9 

V Comparison  of  Results  of  Extensional  and  Inextensional  theories  9 

a)  Buoyant  Shell,  ^ a = 1,  n = 1 9 

2m 

b)  Light  Shell,  LI  = .A,  n = 1 9 

2m 

c)  Light  Shell,  JLL  . S,  n = 2 9 

2m 

VI  Shell  with  an  Additional  Mass  11 

VII  Summation  of  Effects  from  all  Modes,  Pressure  of  Medium  12 

VIII  Total  Acceleration,  Stress  and  Fluid  Pressure  for  Small  Values  of  t 13 

References  17 

Appendix  19 

Fig.  1-23  


0,r 
w,  v 


A 

u 

c 


cs 

<1 


K 

111 

mn,  mn 

M 


ptfl.t) 

r 

iyo 

p u> 

n 

V^n 

T 

i 

Vn“> 

X 

Xn 

Z(0,t) 


« 2=  cos 


6 


LIST  OF  SYMBOLS (1) 

Polar  coordinates,  See  I'’ig.  1. 

Radial  and  tangential  displacements  of  the  cylindrical  shell. 
Note  that  a positive  displacement  w is  inward. 

Cross  sectional  area  of  the  cylindrical  shell. 

Radius  of  the  cylinder. 

Velocity  of  sound  in  medium. 

Velocity  of  sound  in  shell. 

Distance  from  neutral  axis  of  the  shell  to  the  extreme  fiber. 
Coefficient  defining  the  shape  of  the  modes  in  vacuo. 

See  Eq.  (III-4)  - (III-6). 

Modulus  of  Elasticity. 

Thickness  of  shell. 

Hankel  function  of  second  kind  of  order  n. 

Modified  liessel  function  of  the  first  kind  of  order  n. 

Spring  constant  per  unit  area  of  shell  (Section  VIII) 

Mass  of  shell  per  unit  area. 

Generalized  masses. 

Additional  mass  per  unit  length  of  shell  (Section  VIII). 
Number  of  circumferential  waves  of  mode. 

Shock  wave  pressure  on  surface  of  the  shell. 

Constant  peak  pressure  of  shock  wave, 
t dependent  expansion  coefficient  of  p(0,t). 
t dependent  expansion  coefficient  of  / (0, t). 

Generalized  coordinates. 

Kinetic  Energy. 

Time. 

Radial  component  of  particle  velocity  of  the  shock  wave. 

t dependent  expansion  coefficient  of  v(#,t). 

Displacement  of  the  added  mass  M. 

Coefficient,  x = X n e®1 
n ’n 

Fluid  pressure  on  the  surface  of  the  shell. 

^(1  - £i)  Angle  of  envelopment 
Sec  Eq.  (III-l). 

Hoop  strain. 


(1)  Additional  symbols  in 


the  appendix  are  defined  as  they  occur. 
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I Introduction 

'Hie  response  of  an  elastic  cylindrical  shell  to  a transverse  shock  wave  was  obtained  by  Mindlin 
and  Illeich  (1)  by  an  expansion  of  the  response  into  normal  modes  of  oscillation,  neglecting  extensional 
effects,  except  in  the  dilational  mode,  n = 0.  An  alternative  solution  of  the  problem  by  Carrier  (2)  also 
neglects  these  extensional  effects.  The  first  part  of  the  present  paper,  Sections  II  — V,  is  devoted  to  the 
presentation  of  a refined  method  allowing  for  the  extensional  effects  in  all  modes.  Although  these  effects 
damp  out  ultimately,  they  are  necessarily  of  importance  very  shortly  after  the  arrival  of  the  shock  wave. 
While  the  largest  stresses  in  the  shell  occur  at  late  stages  when  the  shell  is  fully  enveloped,  shock  ef- 
fects develop  much  earlier  and  the  refined  theory  is  of  importance  for  their  study. 

Except  for  the  inclusion  of  extensional  effects  the  basic  assumptions  and  method  of  approach  are 
identical  with  those  of  reference  1.  The  cylindrical  shell  is  elastic,  circular  in  cross-section,  homogeneous 
and  infinitely  long.  It  is  submerged  in  an  acoustic  medium,  i.e.  one  whose  motions  conform  to  the  linear 
theory  of  waves  of  expansion.  The  shell  is  subjected  to  a plane  shock  wave,  with  its  wave  front  parallel 
to  the  axis  of  the  cylinder.  By  means  of  an  approximation,  the  potential  of  the  diffraction  and  radiation  is 
eliminated  between  the  simultaneous  equations  which  couple  the  generalized  coordinates  of  the  shell  and 
the  fluid.  As  a result,  pairs  of  second  order  linear  simultaneous  differential  equations  are  obtained,  link- 
ing the  essentially  inextensional  bending  modes  of  the  shell  with  the  essentially  extensional  modes. 

The  solution  of  this  system  of  simultaneous  differential  equations  is  reduced  in  general  to  a sum 
of  four  definite  integrals  during  the  time  of  transit  of  the  shock  wave  across  the  shell,  and  for  a step 
shock  wave,  it  can  be  expressed  in  terms  of  tabulated  functions  thereafter. 

The  differences  in  the  results  of  the  present  paper  and  of  the  method  of  ref.  1 are  evaluated  for 
typical  cases. 

Section  VI  contains  a different  generalization  of  the  problem  treated  in  ref.  1.  It  considers  the 
case  of  an  inextensional  shell  which  encloses  and  supports  an  additional  concentrated  mass  on  the 
centerline  of  the  shell. 

The  final  sections  VII  and  VIII  are  devoted  to  a consideration  of  the  total  effects  from  all  general- 
ized coordinates  (modes)  of  differing  n.  It  is  shown  how  quite  simple  approximate  solutions  for  the  total 
effects  can  be  obtained  for  the  early  stages  of  the  interaction  problem.  It  is  found  that  these  solutions  are 
upper  limits  for  the  correct  ones.  It  is  concluded  that  these  earlv  time  solutions  are  significant  for  accelera- 
tions and  shock  effects,  while  the  stresses  due  to  shock  waves  of  long  duration  can  be  determined  by  the 
mode  approach. 
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II)  Interaction  between  Shell  and  Medium. 

Fhe  response  of  an  elastic  cylindrical  shell  to  a trarsverse  shock  wave  is  studied  in  this  paper  by 
considering  the  shell  without  the  fluid  as  a separate  structure  responding  to  the  dynamic  forces  exerted  by 
the  surrounding  infinite  acoustic  medium.  Using  the  modes  of  free  vibration  of  the  shell  in  vacuo  as  gen- 
eralized coordinates,  its  response  can  be  expressed  in  terms  of  the  infinite  number  of  these  modes. 

The  shock  wave  is  assumed  to  be  plane,  with  its  wave  front  parallel  to  the  axis  of  the  shell.  Such 
a wave  will  only  excite  those  modes  of  the  shell  in  which  the  generators  of  the  cylinder  remain  straight, 
and  the  frequencies  and  modes  required  are  those  of  a thin  circular  ring  determined  in  the  Appendix.  The 
individual  modes  can  be  classified  by  the  integral  number  n of  circumferential  waves  . For  each  number 
n 0 there  exist  two  frequencies,  one  corresponding  to  a primarily  inextensional  (bending)  mode  and  the 
other  to  a primarily  extensional  mode.  The  displacements  w and  v.  Fig.  1,  of  the  former  are  given  by: 


w — cos  n 0 
v = -q-  sin  n B 

K 


(n  4 0) 


(II-l) 


while  those  of  the  latter  are  given  bv: 


w = cos 


n B 


(rf=0,  1,  2,  ..) 


(n-2) 


v = -dnsin  n $ 


The  displacement  of  the  shell  mav  now  be  expressed  bv  generalized  coordinates  q^ , qnand 

qn  (n  # 0): 

w (tf,t)  - qn  4 i (qn  + Ifn)  cos  n B (11-3) 


v (B  t)  = i.  (—  - d q ) sin  n B 

n^!  ,in 

T he  fluid  pressure  / acting  on  the  surface  of  the  cylinder  is  expanded: 


( 11-4) 
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n - 0 ' n 

(t)  cos  n B . 

(11-5) 
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Fig.  1-23 


LIST  OF  SYMBOLS*1* 


Polar  coordinates,  See  Fig.  1. 

Radial  and  tangential  displacements  of  the  cylindrical  shell. 
Note  that  a positive  displacement  w is  inward. 

Cross  sectional  area  of  the  cylindrical  shell. 

Radius  of  the  cylinder. 

Velocity  of  sound  in  medium. 

Velocity  of  sound  in  shell. 

Distance  from  neutral  axis  of  the  shell  to  the  extreme  fiber. 
Coefficient  defining  the  shape  of  the  modes  in  vacuo. 


See  Eq.  (III-4)  - (III-6). 
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mn,  mn 

M 


Modulus  of  Elasticity. 

Thickness  of  shell. 

Hankel  function  of  second  kind  of  order  n. 

Modified  Dessel  function  of  the  first  kind  of  order  n. 
Spring  constant  per  unit  area  of  shell  (Section  Vni) 

Mass  of  shell  per  unit  area. 

Generalized  masses. 

Additional  mass  per  unit  length  of  shell  (Section  VIII). 
n Number  of  circumferential  waves  of  mode. 

p(d,t)  Shock  wave  pressure  on  surface  of  the  shell, 
p Constant  peak  pressure  of  shock  wave. 

t dependent  expansion  coefficient  of  p(#,t). 
t dependent  expansion  coefficient  of  Z (fl,t). 

Generalized  coordinates. 

Kinetic  Energy. 

Time. 

Radial  component  of  particle  velocity  of  the  shock  wave. 

t dependent  expansion  coefficient  of  v(#,t). 

x Displacement  of  the  added  mass  M. 

X Coefficient,  x = X q e^1 

n n^n 

Z(0,t)  Fluid  pressure  on  the  surface  of  the  shell. 

« = cos’ll  - ~)  Angle  of  envelopment 
8,8  See  Eq.  (III-l). 

e Hoop  strain. 
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P (t) 
n 

T 

t 

v(0,t) 

V„(t) 


(1)  Additional  symbols  in  the  appendix  are  defined  as  they  occur. 
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I Introduction 

The  response  of  an  elastic  cylindrical  shell  to  a transverse  shock  wave  was  obtained  by  Mindlin 
and  Bleich  (1)  by  an  expansion  of  the  response  into  normal  modes  of  oscillation,  neglecting  extensional 
effects,  except  in  the  dilational  mode,  n = 0.  An  alternative  solution  of  the  problem  by  Carrier  (2)  also 
neglects  these  extensional  effects.  The  first  part  of  the  present  paper,  Sections  II  — V,  is  devoted  to  the 
presentation  of  a refined  method  allowing  for  the  extensional  effects  in  all  modes.  Although  these  effects 
damp  out  ultimately,  they  are  necessarily  of  importance  very  shortly  after  the  arrival  of  the  shock  wave. 
While  the  largest  stresses  in  the  shell  occur  at  late  stages  when  the  shell  is  fully  enveloped,  shock  ef- 
fects develop  much  earlier  and  the  refined  theory  is  of  importance  for  their  study. 

Except  for  the  inclusion  of  extensional  effects  the  basic  assumptions  and  method  of  approach  are 
identical  with  those  of  reference  1.  The  cylindrical  shell  is  elastic,  circular  in  cross-section,  homogeneous 
and  infinitely  long.  It  is  submerged  in  an  acoustic  medium,  i.e.  one  whose  motions  conform  to  the  linear 
theory  of  waves  of  expansion.  The  shell  is  subjected  to  a plane  shock  wave,  with  its  wave  front  parallel 
to  the  axis  of  the  cylinder.  By  means  of  an  approximation,  the  potential  of  the  diffraction  and  radiation  is 
eliminated  between  the  simultaneous  equations  which  couple  the  generalized  coordinates  of  the  shell  and 
the  fluid.  As  a result,  pairs  of  second  order  linear  simultaneous  differential  equations  are  obtained,  link- 
ing the  essentially  inextensional  bending  modes  of  the  shell  with  the  essentially  extensional  modes. 

The  solution  of  this  system  of  simultaneous  differential  equations  is  reduced  in  general  to  a sum 
of  four  definite  integrals  during  the  time  of  transit  of  the  shock  wave  across  the  shell,  and  for  a step 
shock  wave,  it  can  be  expressed  in  terms  of  tabulated  functions  thereafter. 

The  differences  in  the  results  of  the  present  paper  and  of  the  method  of  ref.  1 are  evaluated  for 
typical  cases. 

Section  VI  contains  a different  generalization  of  the  problem  treated  in  ref.  1.  It  considers  the 
case  of  an  inextensional  shell  which  encloses  and  supports  an  additional  concentrated  mass  on  the 
centerline  of  the  shell. 

The  final  sections  VII  and  VIII  are  devoted  to  a consideration  of  the  total  effects  from  all  general- 
ized coordinates  (modes)  of  differing  n.  It  is  shown  how  quite  simple  approximate  solutions  for  the  total 
effects  can  be  obtained  for  the  early  stages  of  the  interaction  problem.  It  iB  found  that  these  solutions  are 
upper  limits  for  the  correct  ones.  It  is  concluded  that  these  early  time  solutions  are  significant  for  accelera- 
tions and  shock  effects,  while  the  stresses  due  to  shock  waves  of  long  duration  can  be  determined  by  the 
mode  approach. 


II)  Interaction  botwoon  Shod  and  Medium. 

The  response  of  an  elastic  cylindrical  shell  to  a transverse  shock  wave  is  studied  in  this  paper  by 
considering  the  shell  without  the  fluid  as  a separate  structure  responding  to  the  dynamic  forces  exerted  by 
the  surrounding  infinite  acoustic  medium.  Using  the  modes  of  free  vibration  of  the  shell  in  vacuo  as  gen- 
eralized coordinates,  its  response  can  be  expressed  in  terms  of  the  infinite  number  of  these  modes. 

The  shock  wave  is  assumed  to  be  plane,  with  its  wave  front  parallel  to  the  axis  of  the  shell.  Such 
a wave  will  only  excite  those  modes  of  the  shell  in  which  the  generators  of  the  cylinder  remain  straight, 
and  the  frequencies  and  modes  required  are  those  of  a thin  circular  ring  determined  in  the  Appendix.  The 
individual  modes  can  be  classified  by  the  integral  number  n of  circumferential  waves  . For  each  number 
n > 0 there  exist  two  frequencies,  one  corresponding  to  a primarily  inextensional  (bending)  mode  and  the 
other  to  a primarily  extensional  mode.  The  displacements  w and  v.  Fig.  1,  of  the  former  are  given  by: 


while  those  of  the  latter  are  given  by: 


w = cos  n 6 

v = -j-  sin  n 6 
dn 


n 0 


v = -dnsin  n 0 


(n  VO) 


(11-1) 


(it  = 0,  1,  2,  ..) 


(n-2) 


The  displacement  of  the  shell  may  now  be  expressed  by  generalized  coordinates  q , qnand 
q (n  V 0): 

11  OO 

w (0,t)  = qo  + n £ J (qn  + q"n)  cos  n 0 (II-3) 


v (<9,t)  = 


2 (—  - d <f  ) sin  n 0 

n = 1 dn 


(II-4) 


The  fluid  pressure  7.  acting  on  the  surface  of  the  cylinder  is  expanded: 


7 = 2 Q (t)  cos  n 0 . 

n = 0 n 


(11-5) 


The  equations  of  cylindrical  motion  of  the  shell  may  be  written  in  terms  of  the  generalized  coordinates  qn 
and  qn. 


% + <0n  % ' — 
mn 


n = 1,  2,  3 .. 


(11-6) 


£.  +fi T2  q = 

n >n  ~ 

mn 


wh« 


n = 0,  1,  2,  3 


(II-7) 


tin=  m (1  + j^-) 


m„  = m U +dn) 


(II-8) 

(11-9) 
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2 2 

and  <5*  are  the  natural  frequencies  of  the  vibrations  of  the  shell  in  vacuo,  corresponding  to  the  dis- 
placements of  Eq.  (II-l)  and  ( II-2)  respectively.  The  values  of  these  frequencies  are  given  in  the  Appendix, 
Eq.  (A-8). 

The  coefficient  dR  depends  on  the  ratio  of  the  bending  to  the  extensional  stiffness.  For  cases  of 
practical  interest,  and  for  the  lower  modes  (n  —5),  there  is  practically  no  coupling  between  the  bending  and 
extensional  effects  and  the  value  of  dR  is  very  close  to  dR  = n.  In  such  cases,  simplified  expressions  are 
obtained  for  mn,  mn,  ojn  and  2Jn.  (See  Appendix,  Eqs.  A-12,13,23). 

Proceeding  exactly  as  in  ref.  1,  the  pressure  on  the  shell  is  considered  due  to  the  shock  wave, 
and  due  to  the  waves  reflected  and  radiated  from  the  vibrating  shell: 

Z = p{(9,t)  + p j (a,  0,0  (11-10) 

The  velocity  potential  0 of  the  reflected  and  radiated  waves  is  governed  by  the  equation: 

c2V2<£  = 4>  (11-11) 

where  c is  the  velocity  of  sound  in  the  fluid.  With 


p (0, t)  = 2 P (0  cos  n 0 

n=0  n 


(11-12) 


<f>  (r,  0,  t)  = n^Q  0n(r,t)  cos  n 0 

Eq.  (11-10)  yields  the  relation 

^ = Pn  ^ + P (“.*) 

and  Eqs.  (11-6)  and  (II-7)  become,  on  r = a 
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+ O), 


(P„  + P *„) 


n 


(11-13) 


(11-14) 


(n-15) 
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-2  q-  = 
n ^n 


(P„ 


+ P 0n) 


mn 


(11-16) 


The  radial  velocity  of  a point  on  the  shell  must  equal  the  sum  of  the  radial  components  of  the 
particle  velocities  of  the  shock  wave  and  the  diffracted  wave.  Hence: 


Letting 


w 


v (0,0  + 


d <£n  (a, 0,0 
d a 


(11-17) 


v(0,t)  = — L £ Vn  (t)  cos  n 0 
pc  n = 0 n 


(11-18) 


then,  on  r = a 


d tpn  (a,t) 

da 


(n-i9) 


2 


\ 


( 11-20) 


The  potential  function  0n  can  be  eliminated  by  using  the  approximate  expression: 

d t(in  (a,t)  j(9  0n(a,t) 

da  C dt 

derived  in  ref.  1.  The  main  effect  of  the  approximation  is  to  give  the  reflection  and  radiation,  associated 
with  each  mode  of  the  shell,  the  propagation  characteristics  of  plane  waves  rather  than  cylindrical  waves 

Using  Eqs.  (11-14),  (11-19)  and  (11-20),  the  relation 

0n(t)  = pn(t)  + Vn(t)  -pc(4n  + \)  (11-21) 

is  obtained  and  Eqs.  (11-15)  and  (11-16)  become: 


mi  + 


r1  {<*"  + 


9 p„  + v„  s„ 

+ 6>2  qn  = — = — 

n n 


(11-22) 


- Pc/  -v  9 P„  + V s 

qn  + r-  <4n  + Tn)  + 4 % = -=f-  (11-23) 

mn  mn  mn 

The  response  of  the  shell  to  a quite  general  shock  wave  is  obtained  by  integrating  the  set  of 
simultaneous  second  order  linear  differential  equations  with  the  appropriate  Sn  = Pn  + Vn  . The  fluid 
pressure  is  then  obtained  from  Eq.  (11-21).  For  the  case  of  a step  shock  wave,  the  values  of  Pn  and  Vn 
are  given  in  reference  1. 

Ill)  Response  of  Shell  to  a Step  Shock  Wave  for  n ^ 1. 

For  the  case  n = 0,  only  one  equation  is  obtained  which  is  identical  with  the  one  obtained  in 
reference  I.  This  case  is  therefore  not  considered  further. 

For  n > 1,  Eqs.  (II-22,  11-23)  are  integrated  by  the  method  of  variation  of  constants.  In  the  follow- 
ing equations  the  subscripts  n are  omitted  for  simplicity,  except  in  the  terms  n^  and  mn.  With, 

8 =J1C  , H = Pc 

m"  Sn 

Eqs.  (11-22,  11-23)  become: 

q + 8 (d  + q)+  o>2q  = -JL 

m„ 


(III-l) 


(III-2) 


Let 


q + S((|+q)  + ty2q  = -A- 


q = 2 Ci  e ~zkl 
k = 1 K 


T =kli  Ck  »*ke  Zk‘ 


(III-3) 


(111-4) 


3 


( 


\ 


Substitution  into  the  homogeneous  form  of  Eqs.  (III-2,  III-3)  yields  the  characteristic  equation. 


z^  - (5  + 8)  z^  + (o>2  + 52)  z2  - (5  5 2 + S <o2)  z + (<u2  TL?)  = 0 


and  the  value  of  pk  , 


zk  " ^2k  + a 


Hence,  using  Eqs.  (III-5)  and  (III-6),  Eqs.  (III-4)  become: 


q . I C.  e “V 
k = l k 


a zi  - ®2k  + °>L 
q*kS=1  Ci )Cke-2k‘ 


Considering  as  functions  of  t,  and  substituting  Eqs.  (HI-7)  into  Eqs.  (HI-2,  III-3)  the  method 
of  variation  of  constants  leads  to  the  relations: 


ICre-V  = 


| (zk  ~ 5zk  + v2)  ^ 


C]ce-2kt  = 0 


4 . c 

2 - z,  Ci  e~zkt  = — 
1 k K 


t - zk  (z2  - 5zk  + <y2)-Cke  zkl  = 


Eqs.  (in-8)-(III-ll)  are  four  simultaneous  non-homogeneous  algebraic  equations  which  may  be 
solved  for  the  functions  Ck(t).  Integration  of  these  results  gives  the  functions  Ck(t),  which  when  sub- 
stituted into  Eq.  (IH-7)  yield  the  response  of  the  structure.  The  form  of  the  integrals  appearing  in  these 
equations  will  depend  upon  the  nature  of  the  roots  of  Eq.  (IIT-5)  which  may  appear  in  three  combinations: 

a)  four  real  roots 

b)  two  real  and  two  complex  roots 

c)  four  complex  roots. 
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For  small  values  of  n a simplification  may  be  effected.  If  the  frequency  <u  of  the  primarily  inex- 
tensional  mode  is  sufficiently  small,  the  term  in  Eq.  (HI-2)  may  be  neglected.  When  n = 1,  the  inex- 
tensional  frequency  is  already  equal  to  zero.  Hence  Eqs.  (III-2)-(III-3)  become: 


q + 8 (<|  + f)  =JL 


(III-12) 


q + 8 ($  + q)  + <272  q =JL 


(III-13) 


and  the  corresponding  characteristic  equation: 


[z®  - (8  + 8)  z2  + <w2  z - 8w2]  = 0 


separates  into: 


Z j = 0 


(m-i4) 


(ni-15) 


- (8  + 8)  z?  + oj2  zi  - 8&j2  = 0 (k  = 


(m-i6) 


The  corresponding  values  of  are: 


Hence: 


f*k  — 4 - 1 (k  = 2,3,4) 


= ^1  + f Ci  e zkl 


X - 1 (i. 


1)  Cke-Zkl 


(III -17) 


(m-i8) 


The  method  of  variation  of  constants  yields  the  values  of  C^: 


C9  f1  S(X)  e*2*  dA 

z pazn  Jo 


= - J3_  f1 

paz«  Jo 


S(A)  e ^ d\ 


( m-19) 


C4  = - ^4  S(A)eZ4XdA 
pazA  Jo 


_J + _!l  + J!l  r s (a)  dA 

paz2  pazg  pazo  Jo 


1 


where 


[n  mn  1 
ZgS  + z 48  - ZgZ^  - 5^(1  + ) 

J 


mn~®^*2  “ z3>  u4  “ z2> 


j)  (z  A - Zo) 


paz3 


[z2^ 


+ z4S  - z2z4  - 5^(1  + 


i>  1 

J 


»?3  = 


nn  5 ^z4  “ z3^  ^z3  “ z2^ 


*4  “ 


pazi  z25  + z3®  - Z2Z3  - S^(l  + -^) 

L J 


nn  5 (z4  - z3)  (z2  - 24) 


(III-20a) 


(m-20b) 


(III-20c) 


Differentiation  of  Eq.  (Ill- 1 8)  yields  the  velocities  of  the  inextensional  and  extensional  responses: 


i " | - zkcke  Zk‘ 

^ = | - zk  ^ cke_zkt 


(m-21) 


Substitution  of  Eqs.  (III-19)  into  Eqs.  (III-21)  and  a change  of  variables  from  A to: 

y = cos*1  (1  - cA/a) 


(III-22) 


e_z2  f cos  y j y 

e~z3  I cos  7 dy  (III-23) 

e-z4  I c°8  Y dy 

where 

« = cos"1  (1  -— ) (m-24) 

a 

A similar  substitution  and  transformation  in  Eqs.  (111-18)  and  (III-21)  yields  expressions  for  q, "q,  and"q 
respectively. 

2a 

During  the  time  of  transit  (0  £ « £ ;r,  OS  t £ — — ) 0f  the  shock  wave  front  across  the  shell,  the 
integrals  in  Eq.  (III-23)  are  incomplete  Bessel  integrals.  These  functions  do  not  appear  to  be  tabulated. 


yields: 


P c x „z2  -5  cos  <* 

— 1 = 1,'  c 


S(y) 


+ v c 


ZA  a 


3 -S  cos 


+ »?4  « c 


4 ~ cos 


f 

Jo  p 

■ r- 

Jo  p 

c f S . 

J — sin  y 
Jo  p 


sin  y 


sin  y 


6 


\ 


{ 


but  the  integrations  can  be  performed  numerically  without  difficulty.  Upon  and  after  full  envelopment 
(t  ^ the  integrals  have  the  following  forms: 


£y  sinye-ZC08ydy  = i [ez  - !„<*)] 


7*  sin2  y e-ZC09y  dy  =-ll1(z) 


sin  2 X sin  y e z cos  y d y ^ 12  (z) 


T si 

•'IT 


sin  y e—z  cos  y dy  = J_  (e“z  cos  * -ez) 


^ Js„  (y)  .h  y .“I—  V dy  - [l„-  <,t)  - I„  («k)]  „ > 2 


(III-25) 


Hence,  for  t — 2a/c,  the  inextensional  component  of  the  velocity  becomes: 


*£['- 


P c z2a  cos  « 

*1  n = 2—  e -JT 

P z2a 


73c  z3a  cos 
+ 2 e c 


y.  Z i8 

t?4C  *♦  cos  « 


oc  Z ti8  Z f)8  Z 98  "1 

['„<-§->  - 2|i<4>  * '24>iJ 

'[I.  <¥>  - 2'i4>  * '2<?>]] 

z4a  z4a  z4a  q 

l'o  <— > - 2I1<— 1 1 * bMj 
[ {&  - K <¥>  ] 

[zoa  Zoa  1 

,;4)  - u-hJ 

r z4a  24a  1 

> - !n  4>J 


(III-26) 


(III-27) 


Similar  expressions  apply  for  q,  q,  and  q. 

The  various  results  given  are  suitable  for  computations  only  if  Eq.  (III-16)  has  three  real  roots. 
This  equation  may  have  one  real  and  two  complex  roots,  in  which  case  more  suitable  expressions  are  given 
in  reference  3. 
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The  components  of  the  acceleration  are  obtained  from  Eqs.  (III-12)  (III-13)  once  <},  q,  and  q have 
been  evaluated: 


p a pan  p a p c 

2-^n  ' 2~  - " 2—  — (4n  + V 

P ran  P mn  P 


(m-28) 


S - 2 

P a ..  pan  p a pc  9 Pc 

2~  ■ 2FT‘  2T  t ^ * V - ^-r’  — 

p mn  P mn  P c pa 


(III-29) 


where  ■ ■ is  the  average  acceleration  during  the  transit  time,  that  would  produce  a shell  velocity  equal  to 

Z p Q 

the  particle  velocity  in  the  shock  wave. 

The  pressure  of  the  fluid  on  the  shell. 


Z = Q cos  n 6 
n n 


(in-30) 


is  found  by  means  of  Eqs.  (11-21  and  (III-12)  from: 


p p q" 


(m-3i) 


and  is  proportional  to  the  primarily  inextensional  component  of  the  acceleration. 


The  direct  stress  (hoop  stress)  in  the  shell  is  obtained  by  the  relation: 


(1  - v2) 


(III-32) 


1 / d v. 

( ( w + dd ) 


(111-33) 


Using  Eqs.  (II-3)  and  (II-4)  the  n - th  component  of  the  direct  stress  becomes: 


(1  - v )a 


(n  + 1)  qn  cos  n 8 = 


(III-34) 


The  n - th  component  of  the  flexural  stress  is  obtained  from 


an  = ~(n-~  Ed- (q"  + cos  n 6- 


(m-35) 


( 


IV  Accuracy  of  tho  Approximation,  Eq.  (11-20). 

The  approximation,  Eq.  (11-20),  which  haa  been  derived  in  reference  1 is  valid  for  short  times  after 
the  impact  of  the  wave.  It  is  applicable  only  for  such  values  of  the  time  t where 


II  (2)  (J2_)  ~ 5 H /{2)  (— ) . (IV-1) 

n ct  n ct 

The  approximation  is  good  for  large  values  of  the  argument  of  the  Hankel  Functions,  but  it  depends  on  the 
value  of  n how  large  it  must  be. 

For  n = 0 the  approximation  is  still  meaningful  at  the  end  of  the  transit  time.  At  this  time  the 
value  of  the  argument  is  ff/2  and  the  error  is  about  10%.  However,  for  n > 1,  the  error  increases  and  the 
results  at  transit  time  or  later  are  certainly  not  applicable.  This  is  fortunately  not  too  serious  because  the 
major  effects  in  the  higher  modes  occur  within  a rather  short  time. 

To  obtain  a criterion  for  judging  the  accuracy  of  the  numerical  results  for  n = 1 given  in  this 
paper,  one  can  compare  the  values  of  the  approximate  solutions  for  t -»\oo  and  the  asymptotic  values 
obtainable  from  the  solution  of  the  problem  by  integral  transform^)  in  reference  2.  The  theory  of  this 
paper  gives 

,im  = _L_ . (iv -2) 

t « p c 


while  the  value  obtained  by  asymptotic  integration  of  the  results  of  reference  2 is 


(2) 


lim 
t - ® 


P 2 

~P  c 1 + 2m 

pa 


(IV-3) 


The  two  values  agree  only  if  = 1,  i.e.  for  a neutrally  buoyant  shell.  Due  to  the  fact  that  the 

m m 

important  peak  accelerations  occur  at  early  times  it  is  likely  that  the  difference  in  the  velocity  at  t -*  oo 
does  not  affect  these  accelerations. 

v)  Comparison  of  Results  of  Extensional  and  Inoxtonsional  Theories 

To  illustrate  the  extensional  effects,  response  curves  for  n = 1 and  2 were  obtained  and  can  be 
compared  with  the  results  found  from  the  inextensional  theory,  reference  1.  To  be  able  to  combine  the 
response  from  the  various  modes,  it  is  necessary  to  use  consistent  values  of  the  frequencies.  The 
selected  values  representing  steel  shells  in  water,  are 

5 ™ (n2  * I)  — (V-l) 

"3  a 

Two  cases  were  considered  for  n = 1 in  order  to  show  that  the  extensional  effects  are  more 
pronounced  for  light  shells  than  for  buoyant  ones. 


(X)  This  solution  did  not  include  extensional  effects,  but  this  should  be  immaterial. 

(2)  This  result  was  previously  obtained  by  Dr.  W.  T.  Settee  using  quite  simple  impulse  considerations; 

the  asymptotic  integration  which  does  not  present  any  fundamental  difficulty,  was  made  to  confirm  this  result. 
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1,  n - 1 


a)  Buoyant  Shell,  «= 

Zm 

Figure  2 shows  the  extensional  effect  on  the  velocity  component  ^ to  be  small;  it  also  shows 
which  has  no  counterpart  in  the  inextensional  theory.  Fig.  3 compares  the  accelerations  and  the 
fluid  pressures  Qj.  The  combinations  of  the  accelerations  qj  and  qj  give  the  acceleration  of  points  on 
the  surface  of  the  shell: 


w = (cjl  + qj)  cos  0 (V-2) 

For  neutral  buoyancy,  the  fluid  pressure  Q^does  not  change  significantly  due  to  extensional  effects. 
Furthermore,  the  acceleration  cf}  which  describes  the  rigid  body  motion  of  the  entire  shell  is  only  slightly 
affected. 

To  obtain  a comparison  of  the  acceleration  w of  the  surface  of  the  shell,  Fig.  4 gives  w (0,t),  i.e. 
w for  0 = 0.  It  is  seen  that  the  peak  of  the  extensional  response  is  nearly  one  third  larger  than  the 
inextensional  one.  It  should  be  noted  however,  that  the  area  under  the  w curve,  i.e.  the  impulse  imparted 
to  an  element  of  the  shell,  is  about  the  same  for  both  theories. 

Fig.  5 finally  shows  the  ratio  of  the  direct  circumferential  stress  to  the  hoop  stress  o*  = -pa/h 
due  to  a static  pressure  of  intensity  p. 

The  extensional  curves  have  a heavily  damped  oscillatory  character  not  shown  by  those  from  the 
inextensional  theory.  This  is  clearly  due  to  the  influence  of  the  extensional  vibrations  of  the  shell. 

b)  Light  Shell,  = 5,  n = 1 

Figures  6 to  9 show  the  corresponding  curves  for  a shell  whose  mass  is  only  one  fifth  of  that  of 
the  displaced  medium.  It  is  quite  apparent  that  the  differences  between  extensional  and  inextensional 
theory  are  larger  than  that  for  the  buoyant  shell.  All  the  curves  again  have  a heavily  damped  oscillatory 
character  with  the  period  of  the  extensional  vibrations  of  the  shell. 

Figure  6 includes  also  the  asymptotic  value  of  the  velocity  for  large  times  according  to  the 
expression  given  in  Section  IV.  It  is  clear  that  present  theory  breaks  down  at  or  before  the  end  of  the 
transit  time  t = . It  should  be  stressed  that  the  peak  acceleration  and  the  peak  direct  stress  all 

occur  very  much  earlier,  and  they  are  therefore  well  approximated  by  the  present  theory. 

c)  Light  Shell,  ~ = 5,  n = 2 

z m 


Figures  10—14  show  curves  for  the  n = 2 mode  for  the  light  shell  already  considered  in  case  b. 

Fig.  10  shows  the  ratio  of  the  maximum  flexural  stress  o^  to  the  hoop  stress  a*. 

The  extensional  effect  on  the  flexural  stress  is  very  small  and  the  inextensional  theory  curve  is 
indistinguishable  from  that  of  the  present  theory. 

Figure  11  shows  the  extensional  effect  on  the  velocity  component  t0  be  fairly  small,  it  also 
shows  ^2  which  has  no  counteipart  in  the  inextensional  theory.  Figure  12  shows  the  accelerations  and  the 
fluid  pressure  Qg.  The  fluid  pressure  Qg  is  only  slightly  affected  by  the  extensional  effects.  The  acceler- 
ation components  ^ and  ^ ^ *l,em8elve8  have  no  physical  significance  but  their  combination  gives  the 
acceleration  of  points  on  the  surface  of  the  shell: 

w = + qjl  cos  2 0 (V-3) 

Figure  13  shows  the  acceleration  of  the  shell  at  6 = 0,  i.  e.  w (0,t).  The  peak  of  the  extensional  response 
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is  increased  by  about  20%  over  the  inextensional  one.  Again,  the  impulse  imparted  to  an  element  of  the 
shell  is  about  the  same  for  both  theories.  Figure  14  finally  shows  the  ratio  of  the  direct  circumferential 
stress  to  the  hoop  stress  a*. 

Again,  due  to  the  influence  of  the  extensional  vibrations  of  the  shell,  all  the  extensional  curves 
have  a damped  oscillatory  character  not  shown  by  the  inextensional  theory.  The  damping  is  much  less 
pronounced  for  this  mode  than  for  the  mode  n = 1. 

VI.  Shsll  with  an  Additional  Mass 

It  is  possible  to  treat  the  slightly  more  general  case  of  a cylindrical  shell  with  an  additional 
mass  M ( per  unit  length)  on  the  axis  in  a similar  manner  as  the  plain  cylindrical  shell.  The  mass  is 
assumed  to  be  connected  to  the  shell  by  “uniformly  distributed”  springs  having  a spring  constant  K 
(per  unit  of  area  of  the  shell).  Fig.  (15).  To  show  that  the  interaction  problem  leads  to  similar  equations 
as  in  the  case  without  the  mass  M treated  before,  the  modes  in  vacuo  are  required.  For  simplicity  only 
the  case  of  an  inextensional  shell  is  considered.  The  kinetic  and  potential  energies  per  unit  length  of 
shell  are: 


T ♦ <2>  M Mi2 

2 j o 2 

(VI-l) 

V - ^00  + ^0  + ~7T  fo”  COB  0-w)2  d0 

2a 

where  x is  the  displacement  of  the  mass  M.  For  an  inextensional  shell  the  condition 

w = v0  (VI-2) 

must  be  satisfied.  The  differential  equations  of  motion  can  be  obtained  from  Hamilton’s  Principle;  these 
differential  equations  can  be  solved  by  general  expressions 

w = qn  cos  n 9 e'^n* 

v = — q sin  n 0 e'^n*  (VI-3) 

n 

X = e 
n m 

where  it  is  found  that  the  displacement  x vanishes  except  for  n = 1;  in  other  words  the  motion  of  the 
mass  does  not  couple  with  the  modes  for  n / 1 at  all.  The  only  effect  of  the  mass  is  a slight  change 
in  the  frequencies  due  to  the  springs  K.  The  solutions  for  n ^ 2 found  in  reference  1 apply  therefore 
for  the  new  system,  provided  the  frequencies  are  modified;  if  extensional  results  are  included  it  can  be 
shown  similarly  that  the  results  of  Section  III  apply,  except  for  n = 1. 

For  the  case  n = 1 (inextensional)  two  modes  exist,  having  the  respective  frequencies 


(VI-4) 


-2  K 

"l  “T 


, i autt  v 

in'  + "¥ 


i 


The  frequency  - 0 pertains  to  a rigid  body  motion,  similar  to  the  one  found  previously  for  n - 1; 
the  shape  of  the  mode  being 


w m cos  0 , v = sin  0 , Xj  « 1 
The  shape  of  the  mode  for  the  other  frequency  is 

w « cos  0 , v » sin  0 , Xi  =*  - 


(VI-5) 


2a  mn 


(VI-6) 


It  is  a relative  motion  of  the  mass  M and  the  shell,  without  change  of  shape  of  the  shell  The  values 
of  the  generalized  masses  for  these  modes  are: 


ml  - 2m  (1  + 


5j  - 2m  ( 1 + — tt~~) 


(VI-7) 


The  w component  of  the  new  mode  shapes,  Eqs.  (VI-5,6), does  not  differ  from  that  of  Eqs.  (II-l,  2), 
and  the  solution  of  the  interaction  problem  for  the  inextensional  shell  with  the  added  mass  is  formally 
identical  with  the  solution  already  found  for  the  extensional  shell  without  the  mass  in  the  case  n =«  1.  It  is 
simply  necessary  to  use  the  appropriate  values  of  the  frequencies  and  generalized  masses. 

The  extensional  case  could  also  be  treated,  but  leads  to  three  simultaneous  differential  equations 
of  the  type  (II-6,  7).  If  the  frequency  according  to  Eq.  (VI-4)  is  sufficiently  smaller  than  the  extensional 
frequency  of  the  shell  there  will  be  little  coupling,  and  the  acceleration  of  the  mass  M can  be  found  with 
good  approximation  from  the  inextensional  analysis. 


Fig.  16  shows  the  acceleration  X of  the  mass  at  the  center  for  the  case  of  a light  shell, * 5, 
where  the  additional  mass  M = 8rram  is  selected  to  malce  the  structure  neutrally  buoyant;  a rather  stiff 
spring  was  selected  giving  5^  1.00-^-,  about  one  third  of  the  extensional  frequency  <uq  « 3.33-jj- . 

Only  the  portion  of  the  curve  - 1 which  follows  from  the  closed  expressions  (III-26)  is  shown;  the 
portion  < 1 would  have  required  additional  numerical  integrations  which  were  thought  unwarranted, 
as  the  maximum  value  of  x occurs  for-^  > !•  In  spite  of  the  stiff  spring,  it  is  seen  that  the  largest 
acceleration  of  the  mass  M is  considerably  reduced  and  its  occurrence  delayed  compared  to  the  values 
of  the  accelerations  w (0,t)  on  the  surface  of  the  shell;  Fig.  16  indicates  for  this  purpose  the  w curves 
found  previously  for  a light  and  for  a buoyant  shell  without  center  mass  M.  A similar  computation  for 
the  case  m-i  “ 0.10 — shows  a reduction  of  the  acceleration  to  max.  i ~ 0.03  _E_  indicating  the  protec- 

® pa 

tion  against  shock  effects  which  can  be  obtained  by  shock  mountings. 


VII  Summation  of  Effects  from  all  Modes:  Pressure  of  Medium 

To  obtain  the  actual  physical  effect  of  the  shock  wave  it  is  necessary  to  add  the  contributions 
from  all  modes.  While  it  is  correct  that  for  a step  wave  and  for  large  values  of  the  time  t the  lowest  modes 
n - 0 and  1 will  be  dominant,  this  is  not  the  case  for  very  small  values  of  t where  the  higher  modes  re- 
main important.  It  is,  therefore,  of  considerable  interest  to  realize  that  there  is  a possibility  of  actually 
summing  the  contributions  for  different  n.  This  question  will  be  studied  in  the  next  two  sections. 


(1)  This  is  of  course  due  to  the  assumed  inextensionality. 


A very  interesting  result  which  is  obviously  valid  for  any  shape  of  the  structure,  not  only  for 
cylindrical  shells,  can  be  obtained  by  considering  the  pressure  Z of  the  medium  on  the  shell.  Sub- 
stituting Eq.  (11-21)  into  Eq.  (II-5)  gives 

Z - 2 P_(t)  cos  n 0 + S V (t)  cos  n 6 — pc  § (<L  + q_)  cos  n 0 (VII-1) 

n»0  n n=0  “ r n=0  n n 

Each  of  these  sums  can  be  recognized:the  first  one  is  the  pressure  in  the  shock  wave,  p(0,t),  see  Eq. 
(11-12),  the  second  one  is  a multiple  of  the  radial  component  of  the  particle  velocity  in  the  shock  wave, 
pc  v(0,t)  = cos  #p(0,t),  Eq.  (11-18);  the  last  one  is  the  radial  velocity,  W (0,t).  The  fluid  pressure 
is  therefore 


Z (0,t)  = (1  + cos  6)  p (0,t)  -pew  (0,t)  (VO-2) 

This  relation  expresses  the  fluid  pressure  at  a given  point  solely  as  function  of  the  incident  pressure 
p(0,t)  and  of  the  velocity  w (0,t)  of  the  point  at  which  the  pressure  acts.  Within  the  limit  of  validity  of 
this  equation,  i.e.  for  small  values  of  t,  this  relation  permits  the  direct  writing  of  the  differential  equations 
of  motions  of  the  structure;  the  equation  of  motion  of  the  fluid  no  longer  appears. 

The  derivation  of  Eq.  (VI1-2)  is  possible  because  of  the  approximation,  Eqs.  (3-9),  in  reference  1 
which  results  in  the  pressure  Z at  any  point  of  the  shell  becoming  independent  from  whatever  happens  at 
other  points  of  the  shell;  it  can  therefore  not  matter  if  the  structure  at  those  other  points  would  not  be  a 
cylindrical  shell,  and  one  can  use  the  approximation  Eq.  (VII-2)  therefore  on  any  other  submerged  structure. 
Its  use  is  an  enormous  simplification  of  the  task  of  finding  the  response  of  a structure  to  a shock  wave.  The 
result  can  of  course  only  be  valid  for  the  early  part  of  the  interaction  problem. 

When  Eq.  (VII-2)  is  applied  to  problems  of  skew  incidence  of  a shock  wave  on  a cylindrical  shell, 
or  to  evaluate  the  effect  of  a shock  wave  on  the  end  section  of  a shell,  the  angle  0 must  be  clearly  defined. 
It  is  the  angle  between  the  direction  of  propagation  of  the  shock  wave  and  the  inner  normal  of  the  surface. 


VIII  Total  Accoloration,  Stress  and  Fluid  Pressure  for  Small  Values  of  t 

The  expression  (VII-2)  obtained  in  the  previous  section  will  be  used  to  obtain  approximations  for 
the  total  values  due  to  all  modes  of  the  acceleration,  the  direct  stress,  and  the  fluid  pressure  for  small 
values  of  the  time  t. 

Consider  the  motion  of  the  shell,  Fig.  17,  under  the  action  of  normal  forces  Z(0,t)  defined  by  Eq. 
(VII-2).  For  analysis  purposes  a section  of  the  shell  of  unit  length  in  the  direction  of  the  shell  axis  is 
studied.  This  section  forms  a ring,  similar  to  the  one  considered  in  the  Appendix.  For  small  values  of 
the  time  t the  extensional  forces  in  the  ring  will  be  dominant,  while  the  bending  effects  are  very  small. 
As  the  following  analysis  is  only  intended  for  small  values  of  t the  bending  stiffness  will  be  entirely 
neglected  to  simplify  the  computations.  The  equations  of  motion  of  an  elastically  compressible  ring 
without  bending  stiffness  due  to  a system  of  radial  forces  Z(0,t)  are 


EA 
+ -T 


(w  - yp)  - Z 


EA  , 

a 


w = 


(vm-i) 
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Substituting  Eq.  (VII-2)  for  Z the  equations  of  motion  may  be  written  in  the  form 


* + 


P c 
m 


p(fl,t)  (1 


+ cos 


6) 


m 


EA  , 

9 (w  - 

maz 


(vm-2) 


The  last  term  in  the  first  Eq.  (VHI-2)  is  proportional  to  the  hoop  strain  e = -i-  (w  - vg);  it  is 
physically  obvious  that  a short  time  after  the  shock  wave  hits  the  shell  the  strain  must  be  compressive, 

( > 0.  One  can  reason  further  that  the  strain  t will  decrease  as  9 increases,  giving  y/q  - v^  < 0. 

It  follows  from  these  inequalities  that  the  solutions  w and  v of  the  differential  equations  (VIII-2,  3)  with- 
out the  extensional  terms 


- p c z p(#,t)  ( 1 + cos  6) 
m m 

(vm-3) 

7=0 


will  give  upper  limits  for  w,  but  lower  limits  for  v,  or  their  derivatives.  A solution  based  on  Eq.  (VIII-3) 
will  give  therefore  upper  limits  not  only  foi  the  radial  acceleration  w,  but  also  for  the  hoop  strain  and 
stress  which  are  proportional  to  w - v ^ it  will  also  give  an  upper  limit  for  the  pressure  Z (0,t) 
provided  it  is  determined  from  the  first  of  the  Eqs.  (VIII-1). 

It  is  possible  to  refine  the  above  considerations  and  obtain  a criterion  for  the  length  of  time  for 
which  the  approximate  Eqs.  (VIII-3)  give  the  upper  limits  described.  One  significant  fact  typical  for  the 
conditions  for  very  small  values  of  t is  that  the  shell  ahead  of  the  wave  front  A-A,  Fig.  18,  must  be  at 
rest  and  strainless.  This  is  so  because  during  the  time  dt  it  takes  the  wave  front  to  progress  from 
A to  B a compressive  wave  in  the  shell  will  have  progressed  a distance  csdt  from  A towards  B,  where  c8 
is  the  velocity  of  sound  in  the  shell.  If  the  angle  9 is  small  the  compressive  wave  in  the  shell  will  reach 
point  B later  than  the  wave  in  the  fluid  and  the  shell  ahead  of  the  wave  front  will  be  undisturbed.  This 
will  be  the  case  if  the  time  is  so  short  that  0 < 0^  where  the  limiting  angle  is  defined  by 

°L  = 8in_1  <c/cB)  - c/cs  (Vm-4) 

The  limiting  value  tj^  of  the  time  where  this  occurs  is 

*L  =7  d - 008  °l}  - — 2 (2a)  (VIII-5) 

4cg  c 

It  is  now  intended  to  ascertain  that  Eqs.  (VIII-3)  fumishe  the  described  limits  by  showing  that 


(1)  A lower  limit  for  Z (6, t)  is  obtained  if  it  is  determined  from  Eq.  (VII-2).  One  can  also  obtain  a lower  limit  for  w 
by  using  Eq.  (VIII-2)  in  the  form 


- _ p(0,t)  (1  + cos  0)  p c .1  EA  _ _ 

w - c _ £ w (w  - V, 

m m 2 

ma 


e] 
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for  t ^ t.  the  neglected  term  w - \q  is  positive.  Defining  a new  ‘‘local’*  time  i for  each  point  of  the 
shell  such  that  the  arrival  of  the  shock  wave  occurs  at  t ■ 0 


t - t — — (1  - cos  9) 


the  pressure  can  be  defined  by 


p ( 0,t ) 


(o  if  i < 0 

p (x)  if  z — 0 


(vm-6) 


(vm-7) 


where  p(i)  describes  the  shape  of  the  shock  wave. 

Since  for  t < t^  the  points  of  the  shell  are  undisturbed  prior  to  the  arrival  of  the  wave  front,  the 
displacements  w and  v can  be  expanded  for  i > 0 in  a Taylor  series  in  ascending  powers  of  c.  Because  of 
the  rest  conditions  for  r = 0 the  lowest  terms  of  the  series  will  be  c^, 


& J liU 

^ **  g; 
j - 2 J 


(vm-8) 


(vm-9) 


Substituting  Eq.  (VITI-8)  and  Eq.  (VHI-9)  into  Eq.  (VIII-2)  the  functions  f.  (ff)  and  g.  (9)  may  be  eval- 
uated by  equating  coefficients.  Using  terms  up  to  (3,  the  approximate  expression  for  w - \q  becomes 


p(0,t)  (1  + cos  ff)  I t^c  / p a + cs 

2m  1 L aL  '3m 


2 a 4c  9 4«  « 

cos  t)  s sin  ocoso 

T?  + + 


+ cs  sin2fl  t3 


3m  cL/a 


(vm-io) 


Vg  _ 

where  L = 1 sin20 

_2 


For  the  extreme  case  of  9 = 0,  and  t =.  tj^,  Eq.  (VM-10)  becomes 


w _ V.  . pW,t>  '2t  5 1 c2  \ 

9 — V*  r-2~) 


and  w - v$  will  be  positive  provided 


(vm-ii) 


p a < 5cs 


(vni-i2) 


cg*  , 

For  the  case  of  steel  shells  in  water,  -g-  - 10,  the  above  analysis  will  be ’valid  if  t 5 t^  •»  -tt  of 


P a 

transit  time  and  < 25. 

2m 


/ 


Therefore,  for  times  tip  to  about  1/40  of  fall  envelopement  (and  presumably  somewhat  later),  an  upper 
limit  for  the  acceleration  and  stress  can  be  obtained  from  an  analysis  based  on  Cqs.  (VITI-3),  for  all  but 
extremely  light  shells. 

For  the  case  of  a step  wave,  Eq.  (VIII-3)  yields  the  solution 


+ cos  G)  m 

Tc  [#>  C 


-1)  + t 


(vm-i3) 


The  peak  value  of  the  total  acceleration  w is  quite  high  but  w decays  extremely  fast.  For  small  values 
of  0,  6 - 0, 


- 2p  _P_£  t 
" - e nT 
m 


(vm-i4) 


which  is  the  solution  for  the  case  of  a plane  shock  wave  hitting  an  air  backed  flat  plate  of  mass  m. 
The  upper  limit  for  the  direct  (hoop)  stress  in  the  shell  at  0 * 0 is  given  by 


^ 2m  (V)2 

a * p a ~c  pa 


(vm-15) 


The  upper  limits  for  the  acceleration  w and  direct  stress,  although  valid  only  for  small  t,  are  of 
value  in  judging  the  convergence  of  the  mode  approach.  Fig.  19  shows  the  components  of  the  acceleration 
w at  0 = 0 for  n = 0,  1 and  2 for  the  thin  shell  considered  earlier.  While  ii^is  slightly  smaller  than  Wj 
one  can  not  speak  of  convergence,  and  more  terms  are  obviously  required  to  obtain  a meaningful  result. 

This  is  borne  out  by  Fig.  20  which  shows  the  summation  curve  and  the  upper  limit  approximation  for  w 
obtained  from  Eq.  (VITT-14).  The  latter  may  be  used  to  compute  w up  to  about  ct/2a  = 0.05;  for  larger  t 
the  mode  approach  ought  to  give  reasonable  results  without  using  an  excessive  number  of  terms.  It  will  be 
noted  that  the  area  of  the  two  curves  in  Fig.  20,  representing  the  total  impulse,  differs  not  very  much. 

Fig.  21  shows  the  components  of  the  hoop  stress  a for  n = 0,  1 and  2,  the  sum  of  these  com- 
ponents, and  the  upper  limit  for  the  total  stress  for  small  values  of  t.  The  convergence  of  the 
stress  components  is  quite  good  for  ct/2a  > 0.1,  i.e.  for  the  range  where  the  total  stresses 
are  large.  For  smaller  values  of  t the  approximate  upper  limit  is  available,  but  it  is  unimportant  as  the 
largest  stresses  for  a step  shock  of  long  duration  occur  much  later,  when  the  mode  approach  converges 
well.  It  is  quite  obvious  that  the  significant  part  of  the  stress  is  due  to  the  n = 0 and  1 modes. 

The  mode  approach  is  therefore  quite  satisfactory  for  the  stress  problem  for  shock  waves  of  long  duration. 
For  short  duration  shocks,  say  1/20  of  the  transit  time  or  less,  on  the  other  hand  the  approximate  solution 
is  appropriate  as  being  more  accurate  and  much  simpler. 

Fig.  22  shows  the  components  Qn  of  the  fluid  pressure.  Convergence  is  reasonable  except  for 
small  values  of  t.  Fig.  23  shows  the  upper  limit  according  to  Eq.  (VIII-1)  which  indicates  a rapid  decrease 
of  the  pressure  Z «=  2p  for  t = 0.  For  ct/2a  = 0.10  the  upper  limit  and  the  sum  of  the  first  three  com- 
ponents are  about  alike  indicating  that  the  mode  approach  is  reasonable  for  larger  values  oft. 
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APPENDIX 


Frw  Vibrations  of  o Thin  Circular  Ring. 

The  following  derivation  differs  from  the  usual  treatment  by  considering  simultaneously  bending  and 
extensions!  effects.  The  effects  of  shear  and  rotatory  inertia  are  not  considered,  because  they  affect  only 
very  high  modes  not  of  interest  here. 

Let  w be  the  radial  displacement  ( positive  inward)  of  a point  on  the  center  line  of  the  ring  and  v be 
the  tangential  displacement  (positive  in  the  direction  of  increasing  0).  The  kinetic  energy  of  the  vibrating 
system  is: 

T = ^ fn(w2  + *2)ad0  (A-l) 

Z o 

- E 

and  the  potential  energy  is,  (using  E = for  the  plane  strain  solution), 

(1  - v ) 

V =--1—  J2"^  - w)2  d 0 +—  ?n  (vim  + w)^  dO  (A-2) 

Za  o w 2a^  0 

where  I is  the  moment  of  inertia  of  the  cross  section  of  the  ring  with  respect  to  a principal  axis  at  right 
angles  to  the  plane  of  the  ring. 

Using  Hamilton’s  Principle,  the  differential  equations  of  motion  are  obtained  from  the  condition 
that  the  integral 


J = / * (T  - V)  dt 


(A -3) 


shall  be  stationary.  Substituting  Eqs.  (A-l)  and  (A-2)  in  Eq.  (A-3)  and  applying  the  rules  of  the  calculus 
of  variations,  the  following  set  of  simultaneous  partial  differential  equations  is  obtained: 


EA  , . El  . „ . 

maw  + — (v0  - w)  - -j-  (*oeee  + 2ytee  + «)  = o 


EA  , . 

— mav  4 'v  00  - w0'“  0 

a 


Suitable  solutions  of  Eq.  (A-4)  are: 


w = cos  n 0 e*401 


v «=  — sin  n 0 e,CtJt 
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(A-4) 


(A -5) 


\ 


where  n ■>  0,  1,  2,  ....  is  an  integral  number  indicating  the  number  of  circumferential  waves.  Substitution  of 
Eq.  (A-5)  into  Eq.  (A-4)  leads  to  the  two  homogeneous  equations: 


f,i.“  Jy.il>]  3 ”3 

L ma2  ma*  J “ ma 


(A-6) 


nCA  1 

ma2  + d 
ma 


-vl  - » 

ma  J 


Non-vanishing  solutions  of  Eqs.  (A-6)  and  free  vibrations,  exist  only  if  the  determinant  of  Eq.  (A-6)  van- 
ishes. This  leads  to  the  frequency  equation: 


4 2 Tea 

<U  - 0J 

L ma 


<»2  ♦ 1)  3.  (.2  - 1)2 


ma 


♦ ll  H.  . v - D2  - o 


ma^  ma2 


(A-7) 


Eq.  (A-7)  yields  two  frequencies: 


EA  (.2  ♦ 1)  ♦ 


51 

2ma^ 


(n2  - l)2 


]■ 


(A-8) 


? All  („2  (.2  . !)2]2  . B.  EA_  „2(n2  _ lj2 

^ 2ma^ 


2ma^ 


ma^  ma2 


where  o>n  corresponds  to  a primarily  inextensional  mode  and  &>n  to  a primarily  extensions!  mode.  Substitu- 
tion of  these  frequencies  in  Eq.  (A-6)  yields  values  d and  3 corresponding  to  the  two  frequencies.  Because 
of  the  orthogonality  relations  between  the  two  modes  of  frequency  &>n  and  wn,  the  relation  dn  = — j 
exists  where:  . n 


n^  - 1 (1  - nz)‘ 


2\2 


2n 


2n 


2. 


_l_/n2  + 1 2 
+ T/( 


.2\ 4 ,2 


n 


v 2(1  - n2)3  I (1  - nz)  

' + 2 2.  + 2 4*2 

n a A n a A^ 


(A-9) 


The  displacements  w and  v of  the  primarily  inextensional  modes  are  given  by: 


w «=  cos  n 9 

v = — sin  n 9 

dn 


(A-10) 
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while  those  of  the  primarily  extensional  modes  become: 


w = cos  n 9 

v = -1-  sin  n 9 = -dn  sin  n 9. 


(A-ll) 


For  properties  of  practical  interest  and  for  the  lower  modes  (n<  5),  there  is  little  coupling  between 
the  bending  and  extensional  effects  because  a2A  > >1;  Eqs.  (A-8)  and  (A-9)  may  then  be  approximated  very 
closely  by: 


2 El  (1  - n2)  n2 


and 


*4  (n2  + 1) 


5»  + u 


(A-12) 


dn  = n. 


(A-13) 


Determination  of  the  Generalized  Matzet  and  nv  . 

n n 

The  displacement  of  the  shell  is  given  by  Eq.  (II-3)  and  (II-4).  The  kinetic  energy  of  a section 
of  shell  of  unit  width  for  the  n^1  mode  (n  4 0)  is: 


T - 


[<>  +--2>  £ * a * jM„2] 


(A-14) 


and  the  potential  energy  is: 


v ■ ir  ft ‘ 1)2  + + 1)2  " n2’2 <q2  + ?2) 

The  Lagrangian  equations  of  motion  for  small  oscillations  are: 


(A-15) 


j_  (£T.  jy_  F 

dt  n 

where  Fn  is  the  generalized  force  corresponding  to  the  generalized  coordinate  qn. 


(A-16) 


The  force  Fn  can  be  found  by  the  method  of  virtual  work.  Consider  the  n^1  component  of  the  radial 
pressure  on  the  cylinder,  (Eq.  (II-5) 


Zn  = Q„  (t)  cos  n 9 . 

Upon  changing  the  coordinate  qn  by  an  amount  §qn,  the  corresponding  displacement  w is: 

w = 5qn  cos  n 9 


(A-17) 


(A-18) 
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The  work  done  by  the  pressure  Zn  on  the  above  displacement  becomes: 


V. 


f z«  8’» 


cos  n 9 ad 9 


(A- 19) 


which  when  equated  to  Fn  5qn,  yields  the  generalized  force: 


N 


cos  n 9 ad0  ■ I Q„  cos1’  n 9 ad 9 = Qnrra 


P _ I 7 A _ iA  _ I A 2 

0 0 

The  equations  of  motion  are  obtained  by  applying  Eq.  (A-16) 


(1  +T2)e(n  + maff  Q +T2)  “n  qn  = Qn  "a 


(A-20) 


(A-21) 


marr  (1  + d^)  qn  + man  (1  + d^)  5^  = 

9 —9 

where  and  are  the  frequencies  defined  by  Eq.  (A-8).  Eqs.  (A-21)  may  be  written  finally: 

2 0 
% + % % » -=- 
mn 


where 


- _2  - vn 

qn  + <4  % = -=“ 

mn 


mn  = m d + “"Tf) 
n 


mn  = m (!  + dn) 


(A-22) 


(A-23) 
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Figure  1 
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Office  of  the  Chief  of  Engineers 
Assistant  Chief  for  Works 
Department  of  the  Army 
Bldg.  T-7,  Gravelly  Point 
Washington  25,  D.C. 

Attn;  Structural  Branch  (R.L.  Bloor)(l) 

Office  of  the  Chief  of  Engineers 
Asst.  Chief  for  Military  Construction 
Department  of  the  Army 
Bldg.  T-3,  Gravelly  Point 
Washington  25,  D.C. 

Attn:  Structures  Branch 

(M.  F.  Carey)  (1) 

Protective  Construction 
Branch  (I.O.  Thornley)  (1) 
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Army  (cont.) 

U.S.  Army  Waterways  Experiment  Sta, 

P.  0.  Box  631 

Halls  Ferry  Road 

Vicksburg,  Mississippi 

Attn:  Col.  H.  J.  Skidmore  (1) 

Operations  Research  Officer 
Department  of  the  Army 
Ft.  Lesley  J.  McNair 
Washington  25,  D.C. 

Attn:  Howard  Brackney  (1) 

c)  Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn;  OP-31  (1) 

OP-363  (1) 

Director,  David  Taylor  Model  Basin 


Department  of  the  Navy 
Washington  7,  D.C. 

Attn:  Code  720,  Structures 

Division  (1) 

Code  740,  Hi-Speed 
Dynamics  Division  (1) 

Commander 

Portsmouth  Naval  Shipyard 
Portsmouth,  N.  H. 

Attn;  Design  Division  (1) 

Naval  Ordnance  Laboratory 
White  Oak:  Maryland 

RFD  1,  Silver  Spring,  Maryland 
Attn:  Mechanics  Division  (1) 

Explosive  Division  (1) 

Mech.  Evaluation  Div.  (1) 


Commander 

U.  S.  Naval  Ordnance  Test  Station 

Inyokern,  California 

Post  Office  - China  Lake,  Calif. 

Attn:  Scientific  Officer  (1) 

Chief  of  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D C. 

Attn:  Code  P-314  (1) 

CodeC-313  (1) 

Superintendent 

U S.  Naval  Post  Graduate  School 
Annapolis,  Maryland 


Ballistics  Research  Laboratory 
Aberdeen  Proving  Ground 
Aberdeen,  Maryland 
Attn:  Dr.  C.  W.  Lampson  (1) 

Office  of  Chief  of  Ordnance 
Office  of  Ordnance  Research 
Department  of  the  Army 
The  Pentagon  Annex  #2 
Washington  25,  D.C. 

Attn:  ORDTB-PS  (1) 


Chief  of  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.C. 


Director  of  Research 

(2) 

Code  423 

(1) 

Code  442 

(1) 

Code  421 

(1) 

Officer  in  Charge 

Underwater  Explosions  Research  Div. 


Code  290 

Norfolk  Naval  Shipyard 

Portsmouth,  Virginia  (1) 

Director,  Materials  Laboratory 

New  York  Naval  Shipyard 

Brooklyn  1,  New  York  (1) 

Chief  of  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D-C- 
Attn:  Ad-3,  Technical  Library  (1) 

Rec.  P-  H-  Girouard  (1) 

Naval  Ordnance  Test  Station 
Underwater  Ordnance  Division 
Pasadena,  California 
Attn:  Structures  Division  (1) 

Chief  of  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  TD-41,  Technical  Library  (1) 


Officer  in  Charge 

Naval  Civil  Engr.  Research  &. 

Evaluation  Laboratory 
Naval  Station 

Port  Hueneme,  California  (1) 

Commanding  Officer  & Director 
U.  S.  Navy  Electronics  Lab. 

San  Diego  52,  California 

attn:  Dr.  A.  B.  Focke  (1) 
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Commanding  General 
U.  S.  Air  Forces 
The  Pentagon 
Washington  25,  D.  C. 

Attn:  Res.  & Development  Div.  (1) 

Office  of  Air  Research 
Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  Chief  Applied  Mechanics  Group  (1) 

e)  Other  Government  Agencies 

U.  S.  Atomic  Energy  Commission 

Division  of  Research 

Washington,  D.C.  (1) 


Deputy  Chief  of  Staff,  Operations 
Air  Targets  Division 
Headquarters,  U.  S.  Air  Forces 
Washington  25,  D.  C. 

Attn:  AFOIN-T/PV  (1) 


Director,  National  Bur.  of  Standards 
Washington,  D.C. 

Attn:  Dr.  W.  H.  Ramberg  (1) 
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